Potential Autoimmunity Resulting from Molecular Mimicry between SARS-CoV-2 Spike and Human Proteins.

Molecular mimicry between viral antigens and host proteins can produce cross-reacting antibodies leading to autoimmunity. The coronavirus SARS-CoV-2 causes COVID-19, a disease curiously resulting in varied symptoms and outcomes, ranging from asymptomatic to fatal. Autoimmunity due to cross-reacting antibodies resulting from molecular mimicry between viral antigens and host proteins may provide an explanation. Thus, we computationally investigated molecular mimicry between SARS-CoV-2 Spike and known epitopes. We discovered molecular mimicry hotspots in Spike and highlight two examples with tentative high autoimmune potential and implications for understanding COVID-19 complications. We show that a TQLPP motif in Spike and thrombopoietin shares similar antibody binding properties. Antibodies cross-reacting with thrombopoietin may induce thrombocytopenia, a condition observed in COVID-19 patients. Another motif, ELDKY, is shared in multiple human proteins, such as PRKG1 involved in platelet activation and calcium regulation, and tropomyosin, which is linked to cardiac disease. Antibodies cross-reacting with PRKG1 and tropomyosin may cause known COVID-19 complications such as blood-clotting disorders and cardiac disease, respectively. Our findings illuminate COVID-19 pathogenesis and highlight the importance of considering autoimmune potential when developing therapeutic interventions to reduce adverse reactions.

Comparative Analysis of Structural Features in SLiMs from Eukaryotes, Bacteria, and Viruses with Importance for Host-Pathogen Interactions.

Protein-protein interactions drive functions in eukaryotes that can be described by short linear motifs (SLiMs). Conservation of SLiMs help illuminate functional SLiMs in eukaryotic protein families. However, the simplicity of eukaryotic SLiMs makes them appear by chance due to mutational processes not only in eukaryotes but also in pathogenic bacteria and viruses. Further, functional eukaryotic SLiMs are often found in disordered regions. Although proteomes from pathogenic bacteria and viruses have less disorder than eukaryotic proteomes, their proteins can successfully mimic eukaryotic SLiMs and disrupt host cellular function. Identifying important SLiMs in pathogens is difficult but essential for understanding potential host-pathogen interactions. We performed a comparative analysis of structural features for experimentally verified SLiMs from the Eukaryotic Linear Motif (ELM) database across viruses, bacteria, and eukaryotes. Our results revealed that many viral SLiMs and specific motifs found across viruses and eukaryotes, such as some glycosylation motifs, have less disorder. Analyzing the disorder and coil properties of equivalent SLiMs from pathogens and eukaryotes revealed that some motifs are more structured in pathogens than their eukaryotic counterparts and vice versa. These results support a varying mechanism of interaction between pathogens and their eukaryotic hosts for some of the same motifs.

Dynamic, but Not Necessarily Disordered, Human-Virus Interactions Mediated through SLiMs in Viral Proteins.

Most viruses have small genomes that encode proteins needed to perform essential enzymatic functions. Across virus families, primary enzyme functions are under functional constraint; however, secondary functions mediated by exposed protein surfaces that promote interactions with the host proteins may be less constrained. Viruses often form transient interactions with host proteins through conformationally flexible interfaces. Exposed flexible amino acid residues are known to evolve rapidly suggesting that secondary functions may generate diverse interaction potentials between viruses within the same viral family. One mechanism of interaction is viral mimicry through short linear motifs (SLiMs) that act as functional signatures in host proteins. Viral SLiMs display specific patterns of adjacent amino acids that resemble their host SLiMs and may occur by chance numerous times in viral proteins due to mutational and selective processes. Through mimicry of SLiMs in the host cell proteome, viruses can interfere with the protein interaction network of the host and utilize the host-cell machinery to their benefit. The overlap between rapidly evolving protein regions and the location of functionally critical SLiMs suggest that these motifs and their functional potential may be rapidly rewired causing variation in pathogenicity, infectivity, and virulence of related viruses. The following review provides an overview of known viral SLiMs with select examples of their role in the life cycle of a virus, and a discussion of the structural properties of experimentally validated SLiMs highlighting that a large portion of known viral SLiMs are devoid of predicted intrinsic disorder based on the viral SLiMs from the ELM database.

A Phylogenetic Rate Parameter Indicates Different Sequence Divergence Patterns in Orthologs and Paralogs.

Heterotachy-the change in sequence evolutionary rate over time-is a common feature of protein molecular evolution. Decades of studies have shed light on the conditions under which heterotachy occurs, and there is evidence that site-specific evolutionary rate shifts are correlated with changes in protein function. Here, we present a large-scale, computational analysis using thousands of protein sequence alignments from animal and plant proteomes, representing genes related either by orthology (speciation events) or paralogy (gene duplication), to compare sequence divergence patterns in orthologous vs. paralogous sequence alignments. We use sequence-based phylogenetic analyses to infer overall sequence divergence (tree length/number of sequences) and to fit site-specific rates to a discrete gamma distribution with a shape parameter α. This inference method is applied to real protein sequence alignments, as well as alignments simulated under various models of protein sequence evolution. Our simulations indicate that sequence divergence and the α parameter are positively correlated when sequences evolve with heterotachy, meaning that inferred site rate distributions appear more uniform as sequences diverge. Divergence and α are also positively correlated in both orthologous and paralogous genes, but the average increase in α (as a function of divergence) is significantly higher in paralogous protein alignments than in orthologous alignments. This result is consistent with the widely held view that recently duplicated proteins initially evolve under relaxed selective pressure, promoting functional divergence by accumulation of amino acid replacements, and hence experience more evolutionary rate fluctuations than orthologous proteins. We discuss these findings in the context of the ortholog conjecture, a long-standing assumption in molecular evolution, which posits that protein sequences related by orthology tend to be more functionally conserved than paralogous proteins.

An Easy Protocol for Evolutionary Analysis of Intrinsically Disordered Proteins.

We present an easy protocol for evolutionary analysis of proteins, with an emphasis on studying the evolutionary dynamics of disordered regions. Using the p53 protein family as an example, we provide a guide for finding homologous sequences in a database and refining a dataset before constructing the evolutionary context by building a phylogenetic tree. We show how a multiple sequence alignment and phylogeny for a protein family can be further partitioned into smaller datasets in order to investigate the changes in disorder content across the phylogeny. Based on the evolutionary context, we also investigate site-specific conservation of disorder. Last, we address how to evaluate the evolutionary dynamics of disorder-to-order transitions.

Potential RNA-dependent RNA polymerase inhibitors as prospective therapeutics against SARS-CoV-2.

Introduction. The emergence of SARS-CoV-2 has taken humanity off guard. Following an outbreak of SARS-CoV in 2002, and MERS-CoV about 10 years later, SARS-CoV-2 is the third coronavirus in less than 20 years to cross the species barrier and start spreading by human-to-human transmission. It is the most infectious of the three, currently causing the COVID-19 pandemic. No treatment has been approved for COVID-19. We previously proposed targets that can serve as binding sites for antiviral drugs for multiple coronaviruses, and here we set out to find current drugs that can be repurposed as COVID-19 therapeutics.Aim. To identify drugs against COVID-19, we performed an in silico virtual screen with the US Food and Drug Administration (FDA)-approved drugs targeting the RNA-dependent RNA polymerase (RdRP), a critical enzyme for coronavirus replication.Methodology. Initially, no RdRP structure of SARS-CoV-2 was available. We performed basic sequence and structural analysis to determine if RdRP from SARS-CoV was a suitable replacement. We performed molecular dynamics simulations to generate multiple starting conformations that were used for the in silico virtual screen. During this work, a structure of RdRP from SARS-CoV-2 became available and was also included in the in silico virtual screen.Results. The virtual screen identified several drugs predicted to bind in the conserved RNA tunnel of RdRP, where many of the proposed targets were located. Among these candidates, quinupristin is particularly interesting because it is expected to bind across the RNA tunnel, blocking access from both sides and suggesting that it has the potential to arrest viral replication by preventing viral RNA synthesis. Quinupristin is an antibiotic that has been in clinical use for two decades and is known to cause relatively minor side effects.Conclusion. Quinupristin represents a potential anti-SARS-CoV-2 therapeutic. At present, we have no evidence that this drug is effective against SARS-CoV-2 but expect that the biomedical community will expeditiously follow up on our in silico findings.

Exploring Evolutionary Constraints in the Proteomes of Zika, Dengue, and Other Flaviviruses to Find Fitness-Critical Sites.

Dengue virus (DENV) challenges vaccine design due to antibody-dependent enhancement (ADE) and evidence suggests that Zika virus (ZIKV) experiences ADE with DENV and West Nile virus (WNV) antibodies. Thus, multiple flaviviruses must be considered when developing novel therapies against ZIKV. We analyzed 42 flavivirus polyproteins in their evolutionary context to identify motifs conserved in sequence with low real-time and evolutionary conformational flexibility, thought to be fitness-critical sites. We also analyzed evolutionary rate-shifts between clades for insight on vector specificity. For mosquito-borne flaviviruses, two conserved motifs were identified within the RNA-dependent RNA polymerase (RdRP), critical for flavivirus genome replication. Clade-specific motifs were identified for the ZIKV+DENV and WNV clades, many of which were also in RdRP. Six sites in motifs for WNV experienced significant evolutionary rate-shifts, suggesting their importance for functional divergence. Overall, some of these motifs are prime candidates as broadly neutralizing antiviral drug targets across different mosquito-borne flaviviruses.

Large-Scale Analyses of Site-Specific Evolutionary Rates across Eukaryote Proteomes Reveal Confounding Interactions between Intrinsic Disorder, Secondary Structure, and Functional Domains.

Various structural and functional constraints govern the evolution of protein sequences. As a result, the relative rates of amino acid replacement among sites within a protein can vary significantly. Previous large-scale work on Metazoan (Animal) protein sequence alignments indicated that amino acid replacement rates are partially driven by a complex interaction among three factors: intrinsic disorder propensity; secondary structure; and functional domain involvement. Here, we use sequence-based predictors to evaluate the effects of these factors on site-specific sequence evolutionary rates within four eukaryotic lineages: Metazoans; Plants; Saccharomycete Fungi; and Alveolate Protists. Our results show broad, consistent trends across all four Eukaryote groups. In all four lineages, there is a significant increase in amino acid replacement rates when comparing: (i) disordered vs. ordered sites; (ii) random coil sites vs. sites in secondary structures; and (iii) inter-domain linker sites vs. sites in functional domains. Additionally, within Metazoans, Plants, and Saccharomycetes, there is a strong confounding interaction between intrinsic disorder and secondary structure-alignment sites exhibiting both high disorder propensity and involvement in secondary structures have very low average rates of sequence evolution. Analysis of gene ontology (GO) terms revealed that in all four lineages, a high fraction of sequences containing these conserved, disordered-structured sites are involved in nucleic acid binding. We also observe notable differences in the statistical trends of Alveolates, where intrinsically disordered sites are more variable than in other Eukaryotes and the statistical interactions between disorder and other factors are less pronounced.

Interaction of Peptide Aptamers with Prion Protein Central Domain Promotes α-Cleavage of PrPC.

Prion diseases are infectious and fatal neurodegenerative diseases affecting humans and animals. Transmission is possible within and between species with zoonotic potential. Currently, no prophylaxis or treatment exists. Prions are composed of the misfolded isoform PrPSc of the cellular prion protein PrPC. Expression of PrPC is a prerequisite for prion infection, and conformational conversion of PrPC is induced upon its direct interaction with PrPSc. Inhibition of this interaction can abrogate prion propagation, and we have previously established peptide aptamers (PAs) binding to PrPC as new anti-prion compounds. Here, we mapped the interaction site of PA8 in PrP and modeled the complex in silico to design targeted mutations in PA8 which presumably enhance binding properties. Using these PA8 variants, we could improve PA-mediated inhibition of PrPSc replication and de novo infection of neuronal cells. Furthermore, we demonstrate that binding of PA8 and its variants increases PrPC α-cleavage and interferes with its internalization. This gives rise to high levels of the membrane-anchored PrP-C1 fragment, a transdominant negative inhibitor of prion replication. PA8 and its variants interact with PrPC at its central and most highly conserved domain, a region which is crucial for prion conversion and facilitates toxic signaling of Aß oligomers characteristic for Alzheimer's disease. Our strategy allows for the first time to induce α-cleavage, which occurs within this central domain, independent of targeting the responsible protease. Therefore, interaction of PAs with PrPC and enhancement of α-cleavage represent mechanisms that can be beneficial for the treatment of prion and other neurodegenerative diseases.

Evolution of intrinsic disorder in eukaryotic proteins.

Conformational flexibility conferred though regions of intrinsic structural disorder allows proteins to behave as dynamic molecules. While it is well-known that intrinsically disordered regions can undergo disorder-to-order transitions in real-time as part of their function, we also are beginning to learn more about the dynamics of disorder-to-order transitions along evolutionary time-scales. Intrinsically disordered regions endow proteins with functional promiscuity, which is further enhanced by the ability of some of these regions to undergo real-time disorder-to-order transitions. Disorder content affects gene retention after whole genome duplication, but it is not necessarily conserved. Altered patterns of disorder resulting from evolutionary disorder-to-order transitions indicate that disorder evolves to modify function through refining stability, regulation, and interactions. Here, we review the evolution of intrinsically disordered regions in eukaryotic proteins. We discuss the interplay between secondary structure and disorder on evolutionary time-scales, the importance of disorder for eukaryotic proteome expansion and functional divergence, and the evolutionary dynamics of disorder.

Avoiding Regions Symptomatic of Conformational and Functional Flexibility to Identify Antiviral Targets in Current and Future Coronaviruses.

Within the last 15 years, two related coronaviruses (Severe Acute Respiratory Syndrome [SARS]-CoV and Middle East Respiratory Syndrome [MERS]-CoV) expanded their host range to include humans, with increased virulence in their new host. Coronaviruses were recently found to have little intrinsic disorder compared with many other virus families. Because intrinsically disordered regions have been proposed to be important for rewiring interactions between virus and host, we investigated the conservation of intrinsic disorder and secondary structure in coronaviruses in an evolutionary context. We found that regions of intrinsic disorder are rarely conserved among different coronavirus protein families, with the primary exception of the nucleocapsid. Also, secondary structure predictions are only conserved across 50-80% of sites for most protein families, with the implication that 20-50% of sites do not have conserved secondary structure prediction. Furthermore, nonconserved structure sites are significantly less constrained in sequence divergence than either sites conserved in the secondary structure or sites conserved in loop. Avoiding regions symptomatic of conformational flexibility such as disordered sites and sites with nonconserved secondary structure to identify potential broad-specificity antiviral targets, only one sequence motif (five residues or longer) remains from the >10,000 starting sites across all coronaviruses in this study. The identified sequence motif is found within the nonstructural protein (NSP) 12 and constitutes an antiviral target potentially effective against the present day and future coronaviruses. On shorter evolutionary timescales, the SARS and MERS clades have more sequence motifs fulfilling the criteria applied. Interestingly, many motifs map to NSP12 making this a prime target for coronavirus antivirals.

Paralog-Specific Patterns of Structural Disorder and Phosphorylation in the Vertebrate SH3-SH2-Tyrosine Kinase Protein Family.

One of the largest multigene families in Metazoa are the tyrosine kinases (TKs). These are important multifunctional proteins that have evolved as dynamic switches that perform tyrosine phosphorylation and other noncatalytic activities regulated by various allosteric mechanisms. TKs interact with each other and with other molecules, ultimately activating and inhibiting different signaling pathways. TKs are implicated in cancer and almost 30 FDA-approved TK inhibitors are available. However, specific binding is a challenge when targeting an active site that has been conserved in multiple protein paralogs for millions of years. A cassette domain (CD) containing SH3-SH2-Tyrosine Kinase domains reoccurs in vertebrate nonreceptor TKs. Although part of the CD function is shared between TKs, it also presents TK specific features. Here, the evolutionary dynamics of sequence, structure, and phosphorylation across the CD in 17 TK paralogs have been investigated in a large-scale study. We establish that TKs often have ortholog-specific structural disorder and phosphorylation patterns, while secondary structure elements, as expected, are highly conserved. Further, domain-specific differences are at play. Notably, we found the catalytic domain to fluctuate more in certain secondary structure elements than the regulatory domains. By elucidating how different properties evolve after gene duplications and which properties are specifically conserved within orthologs, the mechanistic understanding of protein evolution is enriched and regions supposedly critical for functional divergence across paralogs are highlighted.

The Nuanced Interplay of Intrinsic Disorder and Other Structural Properties Driving Protein Evolution.

Protein evolution often occurs at unequal rates in different sites along an amino acid chain. Site-specific evolutionary rates have been linked to several structural and functional properties of proteins. Previous analyses of this phenomenon have involved relatively small datasets and, in some cases, the interaction among multiple structural factors is not evaluated. Here, we present the results of a large-scale phylogenetic and statistical analysis, testing the effects and interactions of three structural properties on amino acid replacement rates. We used sequence-based computational methods to predict (i) intrinsic disorder propensity, (ii) secondary structure, and (iii) functional domain involvement across millions of amino acid sites in thousands of sequence alignments of metazoan proteins. Our results somewhat corroborate earlier findings that intrinsically disordered sites tend to be more variable than ordered sites, but there is considerable overlap among their rate distributions, and a significant confounding interaction exists between intrinsic disorder and secondary structure. Notably, protein sites that are consistently predicted to be both intrinsically disordered and involved in secondary structures tend to be the most conserved at the amino acid level, suggesting that they are highly constrained and functionally important. In addition, a significant interaction exists between functional domain involvement and secondary structure. These findings suggest that multiple structural drivers of protein evolution should be evaluated simultaneously in order to get a clear picture of their individual effects as well as any confounding interactions among them.

Functional Diversification after Gene Duplication: Paralog Specific Regions of Structural Disorder and Phosphorylation in p53, p63, and p73.

Conformational and functional flexibility promote protein evolvability. High evolvability allows related proteins to functionally diverge and perhaps to neostructuralize. p53 is a multifunctional protein frequently referred to as the Guardian of the Genome-a hub for e.g. incoming and outgoing signals in apoptosis and DNA repair. p53 has been found to be structurally disordered, an extreme form of conformational flexibility. Here, p53, and its paralogs p63 and p73, were studied for further insights into the evolutionary dynamics of structural disorder, secondary structure, and phosphorylation. This study is focused on the post gene duplication phase for the p53 family in vertebrates, but also visits the origin of the protein family and the early domain loss and gain events. Functional divergence, measured by rapid evolutionary dynamics of protein domains, structural properties, and phosphorylation propensity, is inferred across vertebrate p53 proteins, in p63 and p73 from fish, and between the three paralogs. In particular, structurally disordered regions are redistributed among paralogs, but within clades redistribution of structural disorder also appears to be an ongoing process. Despite its deemed importance as the Guardian of the Genome, p53 is indeed a protein with high evolvability as seen not only in rearranged structural disorder, but also in fluctuating domain sequence signatures among lineages.

Engineering adenylate cyclases regulated by near-infrared window light.

Bacteriophytochromes sense light in the near-infrared window, the spectral region where absorption by mammalian tissues is minimal, and their chromophore, biliverdin IXα, is naturally present in animal cells. These properties make bacteriophytochromes particularly attractive for optogenetic applications. However, the lack of understanding of how light-induced conformational changes control output activities has hindered engineering of bacteriophytochrome-based optogenetic tools. Many bacteriophytochromes function as homodimeric enzymes, in which light-induced conformational changes are transferred via α-helical linkers to the rigid output domains. We hypothesized that heterologous output domains requiring homodimerization can be fused to the photosensory modules of bacteriophytochromes to generate light-activated fusions. Here, we tested this hypothesis by engineering adenylate cyclases regulated by light in the near-infrared spectral window using the photosensory module of the Rhodobacter sphaeroides bacteriophytochrome BphG1 and the adenylate cyclase domain from Nostoc sp. CyaB1. We engineered several light-activated fusion proteins that differed from each other by approximately one or two α-helical turns, suggesting that positioning of the output domains in the same phase of the helix is important for light-dependent activity. Extensive mutagenesis of one of these fusions resulted in an adenylate cyclase with a sixfold photodynamic range. Additional mutagenesis produced an enzyme with a more stable photoactivated state. When expressed in cholinergic neurons in Caenorhabditis elegans, the engineered adenylate cyclase affected worm behavior in a light-dependent manner. The insights derived from this study can be applied to the engineering of other homodimeric bacteriophytochromes, which will further expand the optogenetic toolset.

Did the prion protein become vulnerable to misfolding after an evolutionary divide and conquer event?

Despite high sequence identity among mammalian prion proteins (PrPs), mammals have varying rates of susceptibility to prion disease resulting in a so-called species barrier. The species barrier follows no clear pattern, with closely related species or similar sequences being no more likely to infect each other, and remains an unresolved enigma. Variation of the conformationally flexible regions may alter the thermodynamics of the conformational change, commonly referred to as the conformational conversion, which occurs in the pathogenic process of the mammalian prion protein. A conformational ensemble scenario is supported by the species barrier in prion disease and evidence that there are strains of pathogenic prion with different conformations within species. To study how conformational flexibility has evolved in the prion protein, an investigation was undertaken on the evolutionary dynamics of structurally disordered regions in the mammalian prion protein, non-mammalian prion protein that is not vulnerable to prion disease, and remote homologs Doppel and Shadoo. Structural disorder prediction analyzed in an evolutionary context revealed that the occurrence of increased or altered conformational flexibility in mammalian PrPs coincides with key events among PrP, Doppel, and Shadoo. Comparatively rapid evolutionary dynamics of conformational flexibility in the prion protein suggest that the species barrier is not a static phenomenon. A small number of amino acid substitutions can repopulate the conformational ensemble and have a disproportionately large effect on pathogenesis.

Metazoan innovation: from aromatic amino acids to extracellular signaling.

Tyrosine depletion in metazoan proteins was recently explained to be due to the appearance of tyrosine kinases in Metazoa. Here, we present a complementary explanation for the depletion of tyrosine, stating the importance of tyrosine in signaling not only as a phosphorylation target but also as a precursor for catecholamines and hormones. Molecules (dopamine, norepinephrine, and epinephrine, and to a lesser extent serotonin and melatonin) critical to metazoan multicellular signaling are also greatly dependent on a supply of tyrosine. These signaling molecules are synthesized in two highly linked pathways specific to metazoans. In addition, the shikimate pathway that non-metazoans use to synthesize the aromatic amino acids is not present in metazoans. These important pathway changes have occurred between Metazoa and other eukaryotes, causing significant changes to tyrosine metabolism and rendering tyrosine crucial for extracellular signaling. In addition, the evolutionary and functional linkage between these two pathways and the resulting implications for neuropathology are discussed.

Rapid evolutionary dynamics of structural disorder as a potential driving force for biological divergence in flaviviruses.

Protein structure is commonly regarded to be conserved and to dictate function. Most proteins rely on conformational flexibility to some degree. Are regions that convey conformational flexibility conserved over evolutionary time? Can changes in conformational flexibility alter protein function? Here, the evolutionary dynamics of structurally ordered and disordered (flexible) regions are investigated genome-wide in flaviviruses, revealing that the amount and location of structural disorder fluctuates highly among related proteins. Some regions are prone to shift between structured and flexible states. Increased evolutionary dynamics of structural disorder is observed for some lineages but not in others. Lineage-specific transitions of this kind could alter the conformational ensemble accessible to the same protein in different species, causing a functional change, even if the predominant function remains conserved. Thus, rapid evolutionary dynamics of structural disorder is a potential driving force for phenotypic divergence among flaviviruses.

The evolution of protein structures and structural ensembles under functional constraint.

Protein sequence, structure, and function are inherently linked through evolution and population genetics. Our knowledge of protein structure comes from solved structures in the Protein Data Bank (PDB), our knowledge of sequence through sequences found in the NCBI sequence databases (http://www.ncbi.nlm.nih.gov/), and our knowledge of function through a limited set of in-vitro biochemical studies. How these intersect through evolution is described in the first part of the review. In the second part, our understanding of a series of questions is addressed. This includes how sequences evolve within structures, how evolutionary processes enable structural transitions, how the folding process can change through evolution and what the fitness impacts of this might be. Moving beyond static structures, the evolution of protein kinetics (including normal modes) is discussed, as is the evolution of conformational ensembles and structurally disordered proteins. This ties back to a question of the role of neostructuralization and how it relates to selection on sequences for functions. The relationship between metastability, the fitness landscape, sequence divergence, and organismal effective population size is explored. Lastly, a brief discussion of modeling the evolution of sequences of ordered and disordered proteins is entertained.

Evolution of structurally disordered proteins promotes neostructuralization.

Protein structure is generally more conserved than sequence, but for regions that can adopt different structures in different environments, does this hold true? Understanding how structurally disordered regions evolve altered secondary structure element propensities as well as conformational flexibility among paralogs are fundamental questions for our understanding of protein structural evolution. We have investigated the evolutionary dynamics of structural disorder in protein families containing both orthologs and paralogs using phylogenetic tree reconstruction, protein structure disorder prediction, and secondary structure prediction in order to shed light upon these questions. Our results indicate that the extent and location of structurally disordered regions are not universally conserved. As structurally disordered regions often have high conformational flexibility, this is likely to have an effect on how protein structure evolves as spatially altered conformational flexibility can also change the secondary structure propensities for homologous regions in a protein family.